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Abstract 

We present results obtained in the calculation of nuclear ground state prop- 
erties in relativistic Hartree approximation using a Lagrangian whose QCD- 
scaled coupling constants are all natural (dimensionless and of order 1). Our 
model consists of four-, six-, and eight-fermion point couplings (contact inter- 
actions) together with derivative terms representing, respectively, two-, three-, 
and four-body forces and the finite ranges of the corresponding mesonic in- 
teractions. The coupling constants have been determined in a self-consistent 
procedure that solves the model equations for representative nuclei simulta- 
neously in a generalized nonlinear least-squares adjustment algorithm. The 
extracted coupling constants allow us to predict ground state properties of 
a much larger set of even-even nuclei to good accuracy. The fact that the 
extracted coupling constants are all natural leads to the conclusion that QCD 
scaling and chiral symmetry apply to finite nuclei. 

21.10.DR; 21.30.Fe; 21.60.Jz; 24.85.+p 
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I. INTRODUCTION 



Relativistic mean field (RMF) models are quite successful in describing ground state 
properties of finite nuclei and nuclear matter properties. They describe the nucleus as a sys- 
tem of Dirac nucleons that interact in a relativistic covariant manner via mean meson fields 
H^H or via mean nucleon fields fiO| . |TT|j whose explicit forms sometimes derive solely from 



the meson field approaches lfL2|| . The meson fields are of finite range (FR) due to meson 
exchange whereas the nucleon fields are of zero range (contact interactions or point cou- 
plings PC) together with derivative terms that simulate the finite range meson exchanges. 
There are a number of attractive features in the RMF-FR and RMF-PC approaches. These 
include the facts that the combined meson and/or nucleon fields account for the effective 
central potentials that are used in Schrodinger approaches and that the physically correct 
spin-orbit potential occurs naturally with magnitudes comparable to the (empirical) ad hoc 
spin-orbit interactions required in Schrodinger approaches. Equally attractive is the fact 
that for relatively few parameters (~ 10) a vast amount of information is obtained: the 
Dirac single-particle wave functions and corresponding energy eigenvalues, the ground state 
mass, the baryon and charge densities together with their moments, and the properties of 
saturated nuclear matter. Furthermore, these quantities are obtained simultaneously in the 
same self-consistent relativistic Hartree (or Hartree-Fock) calculation. 

In this work we use mean nucleon fields constructed with contact interactions (point 
couplings) to represent the system of interacting Dirac nucleons. We choose this approach 
for the following reasons: (a) possible physical constraints introduced by explicit use of 
the Klein-Gordon approximation to describe mean meson fields, in particular that of the 
(fictitious) sigma meson, are avoided and instead the effects of the various incompletely 
understood and higher order processes are assumed to be lumped into appropriate coupling 
constants and terms of the Lagrangian, as explained in Ref. [fL0~l , (b) the use of point cou- 
plings allows not only (standard) relativistic Hartree calculations to be performed, but also 
relativistic Hartree-Fock calculations [13,13 by use of Fierz relations (up to fourth order 



|), and (c) the use of point couplings, because of their success in the Nambu-Jona-Lasino 
model for the low-momentum domain of QCD |16[ , is perhaps the best way to test for natu- 



ralness of the coupling constants in the seminal Weinberg expansion [O highlighting the role 
of power counting and chiral symmetry in weakening N-body forces. That is, two-nucleon 
forces are stronger than three-nucleon forces, which are stronger than four- nucleon forces, 
resulting in a sequence making nuclear physics tractable. If the dimensionless coupling con- 
stants of the corresponding Lagrangian are of order 1 (natural) then QCD scaling and chiral 
symmetry apply to finite nuclei. Finally, (d), the RMF-PC model allows one to investigate 
its relationship to nonrelativistic point-coupling approaches like the Skyrme-Hartree-Fock 
(SHF) approach and the RMF-FR approach to contrast the importance and roles of the 
different features these models have, as well as to obtain new insights. 

Concerning point (c), the aim of this paper is to determine whether QCD scaling and chi- 
ral symmetry apply to finite nuclei and, by their application, to construct a state-of-the-art 
parameterization of the relativistic mean-field point-coupling Lagrangian. In the following 
we will use the term RMF model for both the version having finite range due to meson 
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exchange, which we call RMF-FR, and the point-coupling (contact interaction) version that 
we denote by RMF-PC. 



Concerning point (d), it is important to note here that one can also view RMF-PC 
as an approach that lies in between the RMF-FR approach and the nonrelativistic Skyrme- 
Hartree-Fock (SHF) approach which is also a well-developed self-consistent mean-field model 
that performs very well (for a review see |n§)- Whereas SHF is based upon density- 
dependent contact interactions with extensions to gradient terms, kinetic terms, and the 
spin-orbit interaction, RMF-FR is based upon a coupled field theory of Dirac nucleons and 
effective meson fields treated at the mean-field level, where density dependence is modeled 
by nonlinear meson self couplings and the role of gradient terms is taken over by the fi- 
nite ranges of the mesons. The kinetic and spin-orbit terms are automatically carried in 
both RMF models || . Thus, a comparison of RMF-PC and SHF addresses the differences 
between in-medium Dirac and Schrodinger nucleons, that is, in kinetic and spin-orbit com- 
ponents, whereas a comparison of RMF-PC and RMF-FR addresses the absence vs. presence 
of finite range and the different treatments of density dependence. Herein we will perform 
these comparisons using precisely the same fitting strategy as in recent SHF and RMF-FR 
adjustments JB?|,f3|,pU|] except that here we will in addition be guided by considerations of 



QCD scaling and chiral symmetry, that is, naturalness. 

We regard the present work with contact interactions as a refined relativistic point cou- 
pling model in comparison to our earlier work [|l^,|ll|,^,^] for the following three reasons. 



First, initial work in determining coupling constants in RMF-PC approaches ||T0| , |i"2f found 
a high correlation among the ground state observables used to determine them, particularly 
the total binding energy and the root-mean-square charge radius. Given this fact and the 
presence of quadratic, cubic, and quartic terms in the various densities appearing in the 
Lagrangian (representing two-, three-, and four-body interactions) results in very delicate 
cancellations among the corresponding many-body forces. This means that determination of 
the coupling constants using a nonlinear least-squares adjustment algorithm with respect to 
the corresponding measured ground-state observables is fraught with difficulty because the 
coupling constants are generally underdetermined. Consequently, the search for the mini- 
mum in the chi-squared hypersurface results in the location of many local minima from which 
erroneous conclusions can be drawn. Herein we address this problem more completely by 
applying two different nonlinear least-squares adjustment algorithms and, finally, developing 
a combined adjustment algorithm that is used to determine our present results. Second, in 
our initial work we considered only spherical even-even closed-shell nuclei or closed-subshell 
nuclei in both proton number Z and neutron number N because, due to explicit omission of 
the pairing interaction, we allowed only orbital occupation probabilities of or 1. Here, we 
introduce orbital occupation probabilities for both protons and neutrons through a standard 
BCS approach in which the proton and neutron pairing strengths are simultaneously deter- 
mined with the coupling constants in the adjustment algorithm. Third, most of our earlier 
work addressed the question of naturalness after the fact, that is, without consideration 
of the complete set of 10 possible Lorentz invariants that may occur (scalar, vector, pseu- 
doscalar, axial vector, tensor, and the same coupled to isospin f) and without consideration 
of the QCD mass-scale ordering of the terms of the Lagrangian. The former consideration 
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is neccessary to properly pursue the question of naturalness while the latter consideration 
leads to, among other things, classification of the (allowed) terms of the Lagrangian accord- 
ing to their relative strengths, which, of course, assists in its construction in the first place. 
We address, and remain cognizant, of both of these considerations in the various approaches 
presented here. 

The paper is structured as follows. The Lagrangian of our relativistic point coupling 
model, together with its variants, is given in Sec. II. Included are the corresponding rel- 
ativistic Hartree equations, expressions for the various densities and potentials appearing, 
and expressions for the calculated observables that are to be used in determining the cou- 
pling constants of the Lagrangian. The approximations that we invoke are also stated here. 
In Sec. Ill we describe the determination of the coupling constants using four different 
least-squares adjustment algorithms with respect to well-measured ground state observables 
and the external constraint of always obtaining reasonable calculated values of the prop- 
erties of saturated nuclear matter. A relatively new approach to the x 2 minimization has 
been developed and we explain how and why. Our results are given in Sec. IV. First, we 
present comparisons of calculation and experiment for nuclei whose measured observables 
were used to determine the coupling constants. Second, we present comparisons of predicted 
and measured observables for nuclei not used in determining the coupling constants. Third, 
we compare our results to those of other RMF approaches. Then we give our final nuclear 
matter predictions and we mention initial results obtained in calculating fission potential 
energy surfaces and properties of superheavy nuclei. We address the role of QCD scaling and 
chiral symmetry in Sec. V where we test our final sets of coupling constants for naturalness 
and present the corresponding evidence obtained that QCD and chiral symmetry apply to 
finite nuclei. Our conclusions and intentions for future work are given in Sec VI. 

II. THE MODEL 
A. The Lagrangian 

The elementary building blocks of the point-coupling vertices are two-fermion terms of 
the general type 

(?KW0 , T e{l,7i} , T G {1,7m, 75, 7s7m, ^} (!) 

with ip the nucleon field, Tj the isospin matrices and T one of the 4x4 Dirac matrices. There 
thus is a total of 10 such building blocks characterized by their transformation character in 
isospin and in spacetime. 

The interactions are then obtained as products of such expressions to a given order. 
The products are coupled, of course, to a total isoscalar-scalar term. By "order" we mean 
the number of such terms in a product, so that a second-order term corresponds to a four- 
fermion coupling, and so on. In second order only the ten elementary currents squared and 
contracted to scalars may contribute, but at higher orders there is a proliferation of terms 
because of the various possible intermediate couplings. 
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In analogy to the nonrelativistic Skyrme-force models, one goes one step beyond zero 
range and complements the point- coupling model by derivative terms in the Lagrangian, as 
e.g. dvipTiTjip . The derivative is understood to act on both %p and tp. The derivative terms 
simulate to some extent the effect of finite range and there may be genuine gradient terms 



from a density functional mapping, as appears e.g. in electronic systems [[23 
In the present work we consider the following four-fermion vertices: 

isoscalar-scalar: (ipip) 2 (= a-field) 

isoscalar- vector: ( , ip r y fM ip) (-07^-0) (= cu-field) 

isovector-scalar: (iprip) ■ (tprip) (= 5-field) 

isovector-vector: (^-ry^) • (ipr^ip) (= p-field) 

and their corresponding gradient couplings d u (...)d v (...). 

These constitute a complete set of second-order scalar and vector currents whose coupling 
strengths in the corresponding Lagrangian we wish to test for naturalness. All of them ex- 
cept for the derivative term for isovector-scalar coupling have appeared in previous RMF-PC 
work |[10| , pJ] , |2T| , p2| , p!2[] , however, the isovector-scalar interaction (5-meson exchange) has been 



found not to improve the description of nuclear ground-state observables in RMF-FR work 
|5|,|7]]. We therefore ask whether the insensitivity of the RMF-FR calculations to the absence 
or presence of delta-meson exchange is due to cancelations, other missing terms, unnatural- 
ness, or a symmetry, and we will investigate this same insensitivity in our RMF-PC work 
here. That is, no RMF-PC calculation has yet included simultaneously the four-fermion cou- 
pling plus corresponding derivative for the isovector-scalar field, so we will do so here. We 
postpone tensor couplings and 3rd and/or 4th order mixed couplings [('0-0) (^7^-0) (-07 M, 0) , 
for example] to our next work which will use the results from this work as the starting point. 
For that work it is important to note that whereas tensor couplings have had little effect 
in RMF-FR calculations || they do have noticable effects in recent RMF-PC calculations 
I2|. Finally, the pseudoscalar channel (tt- meson) is not included here because it does not 



contribute at the Hartree level. 

In this work we begin with a set of higher-order terms that are common to existing 
RMF-FR and RMF-PC studies. These are the familiar nonlinear terms in the scalar cou- 
pling, (ipip) 3 and (^'0) 4 , as well as a nonlinear vector term [(ipj^ip) (tpj^ip)} 2 as used in some 



RMF-FR [24] and RMF-PC |T(| models. Finally, of course, the electromagnetic field and 



the free Lagrangian of the nucleon field must be included. 

Combining all of these terms, we obtain the Lagrangian of the point-coupling model as 
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£ = £fee + £ 4t + C hot + £ dcr + £ em , 

£frCC _ ^([ry Ql* — Uljl/j, 

£hot = _l /3s( ^ ) 3_ l 7s (^)4_ I 7v [(^)(^)] 2 , 

£ em = -eA^[(l - r 3 )/2] 7 ^ - J F^F**" . 

Note that we use the nuclear physics convention for the isospin where the neutron is asso- 
ciated with r 3 = +1 and the proton r 3 = —1. 

As it stands this Lagrangian contains the eleven coupling constants as, ay, ckts, otv, 
Ps, 7s, 7v, ^s, 5y, 5ts, and 5tv- The subscripts indicate the symmetry of the coupling: 
"S" stands for scalar, "V" for vector, and "T" for isovector, while the symbols refer to the 
additional distinctions: a refers to four-fermion terms, 5 to derivative couplings, and (3 and 
7 to third- and fourth order terms, respectively 



The model thus contains one or two free parameters more than analogous RMF-FR 
models. This happens because most RMF-FR models make the tacit assumption that the 
masses in the uj- and p-field can be frozen at the experimental values of the really existing 
mesons. The assumption is justified to the extent that the actual fits to observables are not 
overly sensitive to these masses. In the RMF-PC model, however, experience will still have 
to show whether the derivative-term coefficients can be eliminated in a similar way, so that 
for the present work all parameters are regarded as adjustable. 



B. The Mean-Field and No-Sea Approximations 

Similar to the RMF-FR approach, we consider the RMF-PC approach as an effective 
Lagrangian for nuclear mean-field calculations at the Hartree level without anti-nucleon 
states (no-sea approximation). We thus obtain the mean- field approximation 

i)0 T Tij) — -> w a KO T T<p a (3) 

e a >0 



where the w a are occupation weights to be determined by pairing, see Sec. [HE], the 4> a are 



the Dirac four-spinor single-particle wavefunctions with upper and lower components g and 
/, and the e a are the corresponding single-particle energies. The "no-sea" approximation 
is embodied in the restriction of the summation to positive single-particle energies. All 
interactions in the Lagrangian, Eq. (0), are then expressed in terms of the corresponding 
local densities: 
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isoscalar-scalar: 




= Ha 4 

i. — JUL T 


n{r)<t>n,{r), 


isoscalar-vector: 


Pv(r) 


= Ea4 


a(r)7O0a(r), 


isovector-scalar: 


PTs(r) 


= E a 4 




isovector- vector: 


pTv(r) 




a(r)T 3 7 Q (/) a (r^ 


proton: 


Pc(r) 


= \ {Pv 


(f) -pTv(r)] 



C. Equations of Motion 

Following the mean field approximation, the single-particle wavef unctions <fi a remain as 
the relevant degrees of freedom. Standard variational techniques |25| with respect to 4> Q yield 
the coupled equations of motion for the set {4> a }- We are interested here in the stationary 
solution for which the field equations read 

lo£ a (pa = (i ; y-d + m + V s + Vv7o + V ts t 3 + V r TV r 3 7o + Vb ^70)^0, 



V s = a s p s + (3 s pl + 7spl + ^s^Ps, 
V y = a v pv + 7vPv + <5vAp v , 
Vts = QJtsPts + ^ts^Ptsj 

Vtv = O^TVPTV + ^TV^PTV) 

Vq = eA , AA = -An pc, 
with the local densities as given in Eq. (f§). 



(5) 



D. Relation to Walecka-type RMF models 

The coupling between two nucleon densities is mediated by a finite range propagator 
in the RMF-FR approach. We can expand the propagator into a zero-range coupling plus 
gradient corrections in a standard manner. This gives, for example, in the equation of 
motion for the isoscalar-scalar coupling (see Eq. 2.10b of Ref. 0) 

9 a _ . 9 a 1 9 a a /r \ 

Ps ~ ~2~ Ps + — r Ap s • (6) 



—A + m\ ml m 4 . 

«s S s 

The identification with the corresponding parameters of the RMF-PC is obvious by inspec- 
tion of Eqs. (10-11). The reverse relations read 

ml = a s /S s , gl = -al/5 s . (7) 



Let us have a preview of what parameters to expect. Taking, for example, the set NL-Z2 |2f 
from the RMF-FR we have g a = 10.1369 and m a = 493.150 MeV. This yields an expected 
a s = -4.2252 • 10~ 4 MeV -2 and 5 S = -1.7374 • 10~ 9 MeV" 4 . We will see in Section fTT^ 
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how the optimized RMF-PC parameterizations compare with that. We shall also consider 
the similar relations belonging to the to-, 5-, and p-meson, which are 

(8) 
(9) 
(10) 



ay = 


ml 


5y = 


9l_ 

m i 


a TS = 


mj 


ots — 


9 2 S 


axv = 


m 2 p 


Sty = 


m 4 p 



E. Pairing and the Center-of-Mass Correction 

In order to describe deformed and non-closed-shell nuclei reasonably, pairing has to be 
involved. Since in this paper nuclei close to the drip lines will not be considered, the pairing 
model can be kept simple. For reasons of better comparision, we use precisely the same 



pairing recipe as in former RMF-FR calculations ||27|| . Thus we employ BCS pairing with 
a 5-force and use a smooth cutoff given by a Fermi function in the single-particle energies. 
The occupation amplitudes determined by the gap equation 



' a _ eF 2 + A 2j 



a 



where the Fermi energy ep is to be adjusted such that the correct particle number is obtained. 
The single-particle gaps A a are state dependent and are determined as 

A Q = |d 3 x0;(£)A(£)0 Q (f), (12) 

where A(x) is the pair potential. The occupation weights for the densities (§) are then given 
by w a = v 2 a . The pairing prescription introduces the pairing strength parameters Vp and Vn 
for protons and neutrons, respectively. These are fitted to experimental data simultaneously 
with the coupling constants appearing in the Lagrangian. 

A further important ingredient is the center-of-mass correction and it was shown f2"8 



that the actual recipe for that correction has an influence not only on the predictions for 
light nuclei, but also on the predictions for exotic nuclei. Thus we use here the microscopic 
estimate 

(P 2 ) 

E cm = w> (13) 

where P is the total momentum and M the total mass of the nucleus. Again, we take care 
that the same recipe is used in the RMF-FR forces to which we compare. 



S 



F. Observables 



The determination of the coupling constants included in the model is based upon a 
fit to experimental data, which requires precise numerical comparison of calculated and 
experimentally observed quantities. For mean-field models, the most natural quantities 
are bulk observables such as the binding energy and the moments of the various density 
distributions, while other properties such as the single-particle spectrum cannot be related 
to experiment in a similarly precise way. In the following we discuss the calculation of the 
observables used in our present study. 



1. Binding Energy 



The binding energy for a nucleus with Z protons and N neutrons is computed according 



to 



E B = Zm p + Nm Q - E L - R. 



pair 



Rr 



(14) 



The most crucial part is the energy R^ from the mean-field Lagrangian. It is computed in 
standard manner as Rl = J2a v a e a — J d 3 x(&\£\&) with |$) being the BCS ground state. It 
can be rewritten in terms of the local densities as 



- J d 3 x |a s pg + |avPv + ^"tsPts + §«tvPtv 



+ ¥sP 3 s + ihsPs + i7vPv 



+!<5 s PsAp s + |5 v pvApv + |5tsPtsApts + |5tvPtvAptv 



(15) 



2. Charge Radius 



The starting point for all observables of the charge distribution is the charge form factor 
defined by 



F c h(q) = J d 3 x exp (iq ■ x)p 



ch(X) 



(16) 



It is a function of the momentum transfer q = |q| for spherically symmetric charge distribu- 
tions. Note that the charge density is to be obtained from folding the proton and neutron 
densities with the intrinsic charge (and current) distributions of protons and neutrons. We 
use here the same recipes as previously [0,0. The various bulk properties of the charge 
distributions are deduced in the standard manner p9f 
then as 



We calculate the r.m.s. charge radius 



\ 



r^lhri 

F ch (0) q 



d 2 F ch (q) 
' dq 2 



(17) 
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3. Diffraction Radius 



The diffraction radius is obtained by comparison with a square-well distribution having 
total charge Q: 



'47T 

~3 

Its formfactor is given by 



pl x = Q{^R 3 Y\ r<R, 0, r>R. (18) 



= 3Q ^ (19) 

(ji is the spherical Bessel function). The diffraction radius R is determined such that the 
first root go of the actual formfactor coincides with the first root of the box form factor. This 
is the case when q R is equal to the first root xo of ji(x). We obtain thus the diffraction 
radius by 

x 4.4934095 

Rdms — — • (20) 

<?0 <?0 



4- Surface Thickness 

The height of the first extremum for q > of the charge formfactor can be used to obtain 
information about the characteristic surface thickness a of the charge distribution. It can be 
calculated by comparing the charge formfactor to a distribution that is obtained by folding 
a box-like distribution Pb ox (r) over a Gaussian Pg(r) with width a. The formfactor Fb g (q) 
of the resulting density distribution is simply the product of the single formfactors: 

= F* x (q)F°( q ) = F^(q) exp ( - ±*Y). (21) 

The parameter a is determined in such a way that the original charge formfactor and F bg ' CT (g) 
have the same value at the first extremum q e of F^ ox (q). For R, the diffraction radius as 
obtained in the previous subsection is used. The final result for the surface thickness reads 



1 

a = — 

Qe 

(q e is given by x e /Rd ms , x e being the x- value at the first extremum of ji(x)/x). 



2 In 




(22) 



G. Comments on the Numerical Solution 

The coupled mean-field equations of the RMF-PC and RMF-FR models are solved on 
a grid in coordinate space employing derivatives as matrix multiplications in fourier-space. 
The solution that minimizes the energy of the system is obtained with the damped gradient 



iteration method |30| that has been successfully applied in the solution of such problems. 
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III. DETERMINATION OF COUPLING CONSTANTS 



A. The Task 

The RMF-PC model as presented above contains 11 coupling constants that have to be 
determined in a fit to experimental data. We do this by a least-squares fit i.e. minimization 
of 

2 {or - or pt ) 2 m , 

X = ? (AO,) a ' (23) 

where 0^ xpt are the experimental data and Of 1 denote the calculated values. AOi are the 
assumed errors of the observables which empirically express the demands on the accuracy of 
the model for the respective observables. They are in some cases (for example for binding 
energies) larger than the experimental errors. The observables chosen for the adjustment 
of the coupling constants are discussed below, see section 111 C| . The global minimum of x 2 



should correspond to the optimal set of coupling constants. Finding the global minimum 
is, however, a nontrivial and non-straightforward task. To have a good chance to find the 
optimal set of coupling constants, two different fitting algorithms have been combined. 



B. The Fitting Algorithms 

The fitting of the model parameters is done with a combination of a stochastic method 
and a direct method. The direct method is the Bevington Curved Step |[31|| , which had been 



successfully applied in many previous fits within the RMF model [32] and in the Skyrme 



Hartree-Fock model |l9|,20|. It consists of the Levenberg-Marquardt method |33| with an 



additional trial step in parameter space. The neighbourhoud of the supposed minimum is 
modelled by a parabolic expansion. Close to the minimum this expansion method is used, 
while further away the method switches to sliding down along the gradient in parameter 
space. This method proved to be the most effective one for the given problem and has been 
used for the final minimization. 

This fitting algorithm searches in each step for a lower value of \ 2 . Thus it will walk 
towards the minimum which is closest to the initial configuration in parameter space. This 
often lets it get stuck in a local minimum. One way out is to start the algorithm several 
times with different initial configurations. In this work, a different way out has been chosen 
using Simulated Annealing (SA) J33|. This Monte Carlo algorithm does random steps in pa- 
rameter space. For each step, depending on an externally controlled temperature, the new 
configuration is accepted or refused. At high temperatures, steps that increase x 2 have a 
nonzero probability to be accepted. The more the system is cooled down (we used exponen- 
tial cooling), the smaller the probabilities get for acceptance of a configuration which would 
increase \ 2 . If the cooling is tuned to be sufficiently slow, the routine can escape shallow 
minima and settle down in the global minimum. Additionally, inspired by the techniques of 
Diffusion Monte Carlo |34j , a population of solution vectors (walkers) has been used. Each 



walker, if unsuccessful in terms of x ■> nas a finite probability to die, while on the other hand 
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successful walkers can give birth to additional walkers that start their random walk from 
that parent configuration. 

Fig. |I| demonstrates a sample run (the set of observables for this run is set 1 shown 
below). The maximum number of walkers was set to 4. The two upper figures show the 
population walking through parameter space in the directions as (left) and ay (right). 
The strong correlations among these two parameters are nicely illustrated in correlated 
movements of the walkers. The figures below show the x 2 f° r the different configurations as 
well as the external temperature and the number of walkers. Due to failures of some walkers, 
the population decreases (walkers die after around 1600 iterations ). This information is fed 
back into the algorithm, so that at later times the chance for new walkers to be born increases 
again: after additonal 900 iterations new walkers are born. 

Using that population of walkers (two runs with a maximum number of 4 walkers were 
performed) increases the chance to find the desired minimum of \ 2 ■ Simulated annealing is, 
however, quite slow compared to the other routine and needs considerably more iterations 
than the direct methods. 

To get the best from both classes of fitting algorithms, the following mixed fitting proce- 
dure was used to determine the coupling constants: Two populations of parameter vectors 
were propagated through parameter space with simulated annealing, using the smaller set 
of observables denoted below as set 1. After the system had been cooled down sufficiently, 
these parameter vectors were used as starting points for the Bevington Curved Step proce- 
dure using the larger set 2 of observables. 



C. The Choice of Observables 



Two different sets of observables were employed due to the different aim of the minimiza- 
tion procedures. To explore the parameter space for minima with SA, set 1 was used, which 
is identical to the set of observables used for the determination of the coupling constants 
in Ref. |]10| (see table | for the observables and weights). Since the idea here was to locate 
basins around minima, that set of observables proved to be sufficient to indicate these areas. 



For the final minimization, however, the larger set 2 (see table 0) was used. The pairing 
strenghts were adjusted simultanously, which is important in order to obtain a set of coupling 
constants with predictive power comparable to other mean-field approaches. This set of 
observables had successfully been applied before to fit the parameterization NL-Z2 for the 
RMF model [26] (the only exception being that the pairing strenghts for NL-Z2 used for the 
calculations in the present work were adjusted with a larger set of empirical gaps, see Ref. 



D. The Force PC Fl 



Our first RMF-PC calculation is that corresponding to RMF-FR approaches which treat 
exchange of the a, u, and p mesons. Thus we have three linear terms together with three 
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corresponding derivative terms and three higher order terms. The set of nine coupling 
constants emerging from the fitting procedure with the lowest value of x 2 is called PC-F1 
and is shown in Table (TTT|. 



Note that these coupling constants have correlated errors and uncorrelated errors pi 



The uncorrelated error of a parameter is the allowed variation of that isolated parameter 
(while all other parameters are kept fixed) which enhances x 2 just by the value one. The 
parameters have thus to be given with enough digits that the last digit stays below the un- 
correlated error. This rule is obeyed in the above table. The correlated error of a parameter 
is its allowed change, i.e. within x 2 + l; if & U t ne other parameters are readjusted. Correlated 
and uncorrelated error would be the same if a parameter is completely independent from the 
others. In practice, the correlated errors are much larger than the uncorrelated ones, indi- 
cating strong correlations among the parameters. The largest correlations appear between 
as and ay whose sum happens to provide the largest contribution to the nuclear shell model 
potential. The correlated error of <5tv is quite large and shows that the parameter might 
as well have a positive or zero value. It is quite losely determined by the fitting strategy. 
There is an analogous situation for Skryme forces, where some of the isovector terms posses 
only losely determined parameters. The pairing strengths, on the other hand, have small 
discrepancy between correlated and uncorrelated error. This shows that the fit to pairing is 
basically independent from the fit of the mean field properties. 



E. Quality of the Fit 



The total x 2 f° r PC-F1 is 99.1. Additionally, we consider the \ 2 per point, x\t = 2.11, 
where the number of points is the number of observables taking into account in the fitting 
procedure, which is 47 in our case. The \ 2 P er degree of freedom, Xdi = 2.75, where the 
degrees of freedom are caluclated as the difference between data points and the number of 
free parameters, also measures the quality of the force obtained in the fitting procedure. 
These number need a bit more eludication. To that end, we inspect the (dis-)agreement for 
the various fit observables in detail. This is done in Fig. |2| which demonstrates performance 
of the new RMF-PC force PC-F1 and compares it to the RMF-FR force NL-Z2. One sees 
that the binding energy is described most precisely with an average accuracy of 0.25%. The 
radii are reproduced within about 0.5% percent. The surface thickness comes last. But 
mind that the usage of relative errors punishes this quantity which has a comparatively low 
value of 1 fm. Most actual errors stay within these error bands. There are a few exceptions. 
The energies of 40 Ca and Ni isotopes seem to have trouble and the diffraction radius of 112 Sn 
is a bit large. Comparing the average errors between PC-F1 and NL-Z2, we see slightly 
different trends. NL-Z2 is superior with respect to binding energies and surface thicknesses. 
It does, however, perform less well concerning radii. The total x 2 of NL-Z2 is 132.7 which is 
34% larger than that for PC-F1. The overall performance of the point-coupling thus seems 
to be a bit better, although the difference is not too dramatic. 
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F. Exploring Modifications in the Isovector Channel 

The model Lagrangian (S) contains only a bare minimum of isovector terms. This was 
chosen in close analogy to the RMF-FR. There are many more terms conceivable already 
at the given order of couplings. The problem is that the given obervables all gather around 
the valley of stability and contain only little isovector information. Isovector extensions of 
the model are thus not so well fixed by the data. Nonetheless, it is worth exploring those 
extensions in order to check that one is not missing too much in the above standard model. 

1. Isovector-Scalar Terms 

We now test the linear isovector-scalar term with coupling constant ctxs [see Eqs. (4) 
and (13)]. Table [TV] shows the set of ten optimized coupling constants which we call PC-F2. 
The correlated errors of the isovector coupling constants are much larger than in PC-F1 
(see Table |TTp. The x 2 f°r the extended set given was reduced by only 0.6% compared to 
PC-F1. Thus we find that this extension is not well determined by the present set of data. 
It is interesting to note that the sum of «ts + «tv approximately corresponds to the value 
of cktv in the force PC-F1. This may indicate that the overall isovector strength has a well 
defined value, but the detailed splitting between the two terms is not yet well determined. 

2. Nonlinearities in the Isovector-Vector Terms 

Another obvious extension of the model is the lowest order nonlinear term in the 
isovector-vector density, 

£tv = -~7tv[(#^) • (#VV)] 2 • (24) 

The ten optimized coupling constants which we call PC-F3 are shown in Table [V]. The new 
parameter 7tv is characterized by large uncorrelated and correlated errors, and in addition 
the uncertainties in the parameter #tv have increased compared to the force PC-F1. This 
hints that the experimental observables are unable to pin down the new parameter. The 
overall quality is \ 2 = 98.8 which is only 0.3% better than that of PC-F1. This indicates 
that the extension by a nonlinear isovector term is undetermined at the present stage of the 
fits. 

3. An extended set with 11 coupling constants 

As a last test of possible extensions in the isovector- channel we performed a fit includ- 
ing the four isovector parameters ckts, ^ts, axv, ^tv- The emerging set of eleven coupling 
constants is shown in Table [VT| and is called PC-F4. This set has a small negative coupling 
constant in front of the four-fermion isovector-scalar term leading to a small attraction. 
The sum ckts + «tv leads to a value of ~ 3.3 x 10~ 5 MeV -2 , which is quite close to the 
value obtained for «tv i n the set PC-F1. This observation underlines the statement we 
have already made concerning the force PC-F2, where we saw a similar behavior of the 
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extended isovector strength. Due to the large correlated errors, all isovector parameters 
except «tv are compatible with positive or zero values, showing that the isovector channel 
of this effective Lagrangian is not well determined by the data included in the fit. 



G. Comparison with Walecka-type models 



In Section |11 D| , we estimated expected coupling constants from a gradient expansion of 
the finite ranges in RMF-FR. We compare now the values for the various coupling constants 
with values that we can expect from the finite-range RMF model, choosing the interaction 
NL-Z2 for our comparisons. Table |VII| shows the expected values (except for the isovector- 
scalar channel, since the RMF model with NL-Z2 has no 5 meson) together with the values 
taken from NL-Z2. 

Good agreement can be seen for the coupling constants mainly responsible for the nu- 
clear potential, namely, as and ay, which are very similar in each of the RMF-PC forces 
and are somewhat lower than the corresponding RMF-FR values. 



By looking at the results for the corresponding coupling constants 5s and Sy, we realize 
that there are dramatic discrepancies. In none of the interactions does the sign of 5y agree 
with the RMF-FR value. Only in PC-F4 do all signs of the four isovector coupling con- 
stants comply with the expectations from the estimates [Eqs. (19) and (20)]. One has to 
keep in mind, however, that these coupling constants, due to their large correlated errors, 
are not incompatible with zero. The values for axv agree well with the expected value from 
NL-Z2, reflecting about the same asymmetry energy that all RMF-FR forces deliver (see 
the discussion about nuclear matter properties in the next section). 

One may be suspicious that the different mapping of non-linearities spoils the compari- 
son. To countercheck, we performed one more fit including the 11 coupling constants, but 
setting 7v to zero in order to address the different signs of ay and 5y which appear in all 
sets of coupling constants studied. The resulting set of coupling constants still has the same 
signs, which shows that the negative value of 5y is not related to having nonlinearities in 
the isoscalar-vector channel of the effective Lagrangian. We thus are led to the conclusion, 
that the gradient terms in the RMF-PC embody obviously more than just a compensation 
for the finite range. This may indicate that the present RMF-PC Lagrangian is incomplete. 

Altogether, all isovector extensions turned out to improve the fits only very little. Even 
a detailed analysis of the trends along isotopic chains did not show any significant improve- 
ment. Thus, we did not consider additional forces in our present study because they do not 
appear to be well determined with existing observables. Additionally, the \ 2 P er degree of 
freedom is larger for the extended sets compared to PC-F1, showing that at the present 
stage the extended forces do not incorporate real physical improvements. This may change 
for larger sets of observables which include dedicated isovector data. The large uncertainties 
in the isovector coupling constants in the three extended models shows that there is indeed 
sufficient freedom to accomodate new isovector observables. 



15 



IV. RESULTS 



A. Comparisons 

We now check the predictive power of the newly fitted force PC-F1. We do this by looking 
at the performance for a variety of test cases and observables which were not included in 
the fit. We compare the model both to experimental data and to three other relativistic 
mean-field approaches, namely the older point-coupling model PC-LA [ID|, and the two sets 



NL3 H as well as NL-Z2 J§] from the family of RMF-FR models. NL-Z2 had been fitted 



with precisely the same set of data as PC-F1. PC-LA employed a smaller set of data as 
discussed above. NL3 was fitted to binding energies, charge radii and neutron rms radii 
of the nuclei 16 0, 40 < 48 Ca, 58 Ni, 90 Zr, H6,i24,i32g n and 208p b Additionally, nuclear matter 
properties entered into the fit (E/A = -16 MeV, p = 0.153 fm~ 3 , K = 250 MeV, a sym = 
33 MeV). NL-Z2 and NL3 are two state-of-the-art mean-field forces that have been tested 
in a variety of applications. So this selection of forces will give us a well-balanced picture 
of the quality of modern relativistic mean-field forces. In some cases we compare also with 
state-of-the-art Skryme forces, namely the forces SLy6 |37| and SkI3 GO]. SLy6 aims at 



describing extremely neutron-rich systems up to neutron stars together with normal nuclear 
matter and nuclei. SkI3 has a spin-orbit force that in its isovector properties is analogous 
to the nonrelativistic limit of the RMF-FR model and was fitted using the strategy of Ref. 
p0| which is much similar to the strategy and input data used here. 



B. Nuclear Matter 



Table |VH1] shows the bulk properties of symmetric nuclear matter as predicted by the 
different forces. Like the other RMF approaches, PC-F1 has a rather low saturation density 
of around po — 0.15 fm -3 while the Skyrme forces produce the larger po — 0.16 fm~ 3 (which 
is close to the commonly accepted value [|38|]). Additionally, all RMF forces favor a larger 
binding energy at the saturation point. These are systematic differences between the two 
approaches apparent for both types of RMF as compared to SHF. This indicates that these 
trends are not due to a finite-range in RMF-FR but must have other reasons related to 
relativistic kinematics. 



The incompressibility K of PC-F1 is comparable to that of PC-LA and NL3, whereas 
NL-Z2 produces a much smaller value. The larger value of 270 MeV is much closer to the 
commonly accepted 240 MeV while the value of NL-Z2 is far too small. It is interesting to 
note that the large value of K was aimed at in the fit of NL3 while it just emerged from 
the fit for PC-F1. It is also to be remarked that NL3 achieves this large K at the price 
of producing somewhat too small surface thickness. PC-F1, on the other hand, describes 
surface thickness as well as NL-Z2 (see Figure 0) and has a much larger K than NL-Z2. We 
see here a clear difference of the point-coupling versus finite range. This is corroborated by 
the fact that the SHF models are also point-coupling models and do also tend to predict 
incompressibilities in the range of 250 MeV. 
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The symmetry energy a sym has the same large value in all RMF models while SHF results 
stay closer to the commonly accepted values (~ 30MeV). This is a systematic discrepancy 
between RMF and SHF. It is most probably connected to the rather rigid parameterization 
of the isovector channel in RMF. 

The effective mass is consistently small in all RMF models while SHF can cover a broad 
range of values up to m* jm = 1 and even a bit larger, see e.g. [39j. The value of m*/m in 
the RMF is directly related to the strength of the vector and scalar fields which, in turn, 
determines the spin-orbit splitting. There is thus little freedom to tamper with the effective 
mass in RMF unless one alternative means to tune the spin-orbit force (as e.g. a tensor 
coupling) . 



Fig. [3] shows several features of symmetric nuclear matter as function of density p. 
The results are about similar for NL-Z2, NL3, and PC-F1 while PC-LA shows dramatic 
deviations, particularly for p > 0.17 fm~ 3 . The net potential V = Vs — Vy and the effective 
mass m* play a crucial role to determine the spectra of finite nuclei. Thus we have to 
expect somewhat unusual spectral features for PC-LA. At second glance, we see also slight 
differences between the other parameterizations coming up slowly at larger densities. The 
equation of state E/A is less rigid for PC-F1 (correlated with a slightly smaller potential 
V and less suppressed m*). This is consequence of the fact that the density dependence is 
parameterized differently in point coupling and finite range models. 



C. Neutron Matter 



Neutron matter is a most critical probe for the isovector features. It has been exploited 
extensively in the adjustment of the SHF forces |3"7 |. There are, of course, no direct mea- 
surements. But neutron matter is well accessible to microscopic many-body theory such 
that there exist several reliable calculations of its properties. Figure ^ shows the equation 
of state for the four RMF forces and SLy6. The crosses correspond to data from [|[J. We 
confine the comparison to low densities which are relevant for nuclear structure physics. It 
is obvious that all RMF models show a similar trend which, however, differs significantly 
from the "data" and from SLy6. This is a systematic discrepancy which, again, is related 
to the rather sparse parameterization in the isovector channel. 



D. Binding Energies 

1. Isotopic and Isotonic Chains 

In Figures [| and |6], we show the systematics of relative errors on binding energies along 
isotopic and isotonic chains for the two RMF-PC forces and the RMF-FR forces discussed 
here. All nuclei in these figures are computed as being spherical. Note that the scales are 
different for each figure. As guideline we indicate by horizontal dotted lines the average 
error of the models for this observable. 
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Larger errors show up sometimes for small nuclei in the isotopic chains, see Figure |5]. 
The case 40 Ca is notoriously difficult for PC-F1 and small Ni isotopes are a problem for all 
RMF models. The underbinding of 40 Ca may be excused by a missing Wigner energy EL 



But 56 Ni is already overbound and a Wigner energy would worsen the situation. The rea- 
sons for the deviation have to be searched somewhere else, probably it is again an isovector 
mismatch. 



The heavier systems perform much better. They are described within an error of about 
0.4%, with few exceptions. We also see that NL-Z2 performs best in most cases. Some 
slopes and kinks are also apparent in these plots for all forces. They indicate yet unresolved 
isotopic and isotonic trends. Another interesting observation can be made: the structure of 
the curves is, with differences in detail, similar for NL-Z2 and PC-F1 in almost all cases 
(this is most striking for the Sn isotopes). It shows that the fitting strategy (i.e. the choice 
of nuclei and observables) has direct consequences for the trends of the errors. 

A well visible feature are the kinks of the errors which appear at magic shell closures. 
These kinks indicate that the jump in separation energies at the shell closure is too large 
(typically by about 1-2 MeV). This, in turn, means that the magic shell gap is generally a 
bit to large. Some SHF forces solve that problem by using effective mass m*/m = 1. This 
option does not exist in RMF as we have seen above. But there are other mechanisms active 
around shell closures. The strength and form of the pairing can have an influence on the 
kink (= shell gap). Moreover, ground state correlations will also act to reduce the shell gap 
of the mere mean field description. This is an open point for future studies. 

Figure ^ show the relative errors of binding along isotonic chains, assuming again all 
spherical nuclei. Again, there are larger fluctuations for the small nuclei N = 20 and 
N = 28 while the heavier N = 50 and N = 82 stay nicely within the error bounds. But 
the heaviest iV = 126 chain grows again out of bounds at its upper end. Isotonic chains 
are a sensitive test of the balance between the Coulomb field and the isovector channel of 
the effective Lagrangian. All effective forces discussed here produce larger errors compared 
to the experimental isotonic chains, which shows the need for further investigations of this 
property of the RMF models. 



2. Superheavy Elements 

The upper panel of Figure [7| shows the relative errors in binding energies for the heaviest 
even-even nuclei with known experimental masses (compare with a similar figure in Ref. 
||42|| ) . The lower panel delivers as complementing information the ground-state deforma- 
tions expressed in terms of the dimensionless quadrupole momentum /3 2 . The calculations 
were performed with allowing axially symmetric deformation assuming reflection-symmetric 
shapes. The agreement is remarkable. All forces (with some exceptions for PC-LA) produce 
only small deviations which stay well within the given error band. This is a gratifying sur- 
prise because we are here 40-50 mass units above the largest nucleus included in the fit. It 
is to be noted that most SHF forces do not perform so well and have general tendency to 
underbinding for superheavy nuclei ||42||. There are also (small but) systematic differences 
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between the RMF models. NL3 generally overbinds a little while NL-Z2 and PC-F1 tend 
to underbind. All forces show yet unresolved isovector trends. The increase of the binding 
energy with increasing neutron number is too small. These trends were already apparent 
for known nuclei (see the discussion above). The reasons for all these trends are not yet 
understood. Finally, mind the kinks visible for the Z = 98 and Z = 100 isotopes at neutron 
number N = 152 which hint at a small (deformed) shell closure there. 

All forces predict strong prolate ground-state deformations for these superheavy nuclei 
(/3 2 ~ 0.26 — 0.31). The trends look similar for all forces. The largest deformations appear at 
N = 148 and/or N = 150. But there are systematic differences in detail: NL-Z2 has always 
larger ground-state deformations than the other forces, while PC-F1, PC-LA, and NL3 show 
much similar deformations. The difference is probably related to the surface energy: NL-Z2 
has a lower surface energy than NL3. The symbol with error bars at Z/N=102/152 in Figure 
[7| corresponds to the measured ground-state deformation of 254 No |)4*3"||P|| . This deformation 
is overestimated by all forces, PC-LA and NL3 stay within the error bars, though. The error 
ranges from 6 to 13% which is still acceptable. 



E. Fission Barrier of 240 Pu 

Fig. I] shows the fission barrier of 240 Pu computed in axial symmetry allowing for re- 
flection asymmetric shapes (for a discussion of the numerical methods, see reference [ff5|j). 
The experimental values for ground state deformation, barrier, and isomer energy are taken 
from Ref. All forces predict the same ground-state deformation in agreement with 



the experimental value and they all show the typical double humped structure of the fission 
barrier. Also the first barrier (which corresponds to reflection-symmetric shapes) is very 
similar but too large as compared to experiment. That may be a defect of symmetry re- 
strictions. Triaxial degrees of freedom can decrease the calculated barrier by about 2 MeV 
||45|| , which would bring the curves closer the the experimental value. Moreover, the (yet to 
be calculated) zero-point energy corrections will also lower the barriers somewhat |50| . 



Larger differences develop towards the second minimum and further out (where also the 
asymmetric shapes take over). This can be related to the surface properties of the different 
forces. Forces with a high surface energy place the isomeric state higher up than forces 
with lower surface energy. All forces, however, underestimate the experimental value for the 
energy difference of the ground-state and the isomeric state, which is 2.3 MeV. Vibrational 
zero-point energies may still help in case of NL3. But the minima for the other three forces 
are to deep that those small corrections could bridge the gap. 



F. Observables of the Density 

1. Charge Radius, Diffraction Radius and Surface Thickness 

In this section we take a look at the observables which are related to the nuclear charge 
distribution, the r.m.s. and diffraction radii as well as the surface thickness (see Sec. [Ill C| ). 
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In Figure ^ we show results for Sn and Pb isotopes. The experimental data are taken from 
Ref. IESU51B3. 



The r.m.s. and diffraction radii are described generally good. PC-LA yields to large 
diffraction radii in Sn isotopes. NL-Z2 produces a bit to large radii in Pb isotopes. But 
note that all forces reproduce the trends of the r.m.s. radii in lead with its pronounced 
kink at the magic N = 126. It is a known feature RMF-FR models perform very well in 
that respect [^3|,^0j and we see here that the point-coupling models maintain this desirable 
feature. Larger discrepancies are observed for the surface thickness (lowest panel in Fig. 
|9p. All forces have a tendency to underestimate the surface thickness. That is a common 
feature of the RMF models. NL-Z2 and PC-F1 included that observable in the fit and it is 
then no surprise that they yield an acceptable agreement with data. The two other forces 
produce to small surface thickness. The deviation ranges up to 10%. That is outside the 
range which could be explained by possible ground state correlation effects. 



2. Density Profiles and Formfactors 

Fig. |IT] shows the baryon densities, p\{r) in Eq. (|j), for the nuclei 48 Ca and 100 Sn. They 
all display the typical pattern of a box-like distribution with smoothened surface and oscil- 
lations on top ||29|| . The oscillations are an unavoidable consequence of shell structure. 100 Sn 
shows in addition the suppression of center density due to the repulsive Coulomb force. All 
forces produce about the same bulk properties, i.e. the overall extension, center density, and 
surface profile. But there are sizeable differences for the amplitude of the shell oscillations. 
RMF-FR produces more than factor two larger oscillations than RMF-PC (and even that 
is still a bit larger than the experimentally observed oscillations). The reason is that the 
finite-range folding is more forgiving what these oscillations is concerned. It seems that the 
final nuclear potential is fixed by the data to have in all cases about the same profile with 
not too large oscillations. The amplitude of oscillations in the density carries fully through 
to the potentials in case of point-coupling. Thus the model needs to curb down the initial 
amplitude. In finite-range models, however, the densities are smoothened by folding with the 
meson propagator which gives more leeway for oscillations on the density. Comparison with 
experimental oscillations could help to decide between finite-range and zero-range models. 
But just this observable of shell oscillations is heavily modified by all sorts of ground state 
correlations |53fl. These have first to be fully understood before drawing conclusions on the 



range of the effective Lagrangian. 

For the nucleus 48 Ca, for which experimental data are available, we compare the charge 
formfactor with the predictions of our models. The experimental data are taken from Ref. 



55| , where the charge density is parameterized by a Fourier-Bessel series with the coefficients 
determined directly from the data. This density is then Fourier transformed to obtain the 
formfactor. We show it in Fig. 11, together with the RMF predictions, in the momentum 
range covered by the original analysis. Of special importance are the first root and the height 
of the first maximum for finite momentum transfer, as they correspond to the diffraction 
radius and the surface thickness. We see that all forces overestimate somewhat the first root 
of the formfactor leading to a slightly too small diffraction radius. They reproduce well the 
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following minimum which leads to an accurate prediction of the surface thickness. Note, 
however, that both observables were part of the fitting procedure for the forces PC-F1 and 
NL-Z2. Going to higher momentum transfer, we see that all forces reproduce the second zero 
of the formfactor and that the two RMF-PC forces agree nicely with experiment concerning 
the following maximum, while the two RMF-FR forces overestimate it somewhat. This 
indicates that the momentum expansion of the RMF-PC model appears to work well in that 
respect up to momentun-transfer q « 3.0 [fm -1 ]. 



G. Spin-Orbit Splittings 



Figure 12 shows the relative errors for a selection of spin-orbit splittings in 16 O, 132 Sn, 
and 208 Pb. We have taken care to choose splittings which can be deduced relieably from 
spectra of neighbouring odd nuclei |j56| . All RMF forces, except for PC-LA, perform very 
well. It was shown in a former study that RMF-FR is much better in that respect than 
many Skyrme forces |p6| . We see now that the well fitted point coupling model PC-F1 does 



as well as RMF-FR. The ability to describe the spin-orbit force correctly is thus a feature 
of the relativistic approach. 

The force PC-LA falls clearly below the others. The poor performance is related to 
the too weak fields at large densities, see Fig. |3| and related discussion. The example 
demonstrates that one needs a sufficiently large set of observables to pin down the nuclear 



mean field sufficiently well. The argument is corroborated by Fig. [T3| where we have a quick 
glance at the effective spin-orbit potentials oc V[2mjv — Vs — Vy] _1 . 

The three well performing models have all very similar potentials whereas PC-LA has a 
10% stronger spin-orbit potential which is shifted a little bit to larger radii. This difference 
yields the observed mismatch in the spin-orbit splittings. In turn, this figure shows that the 
allowed variations on the mean fields (here the spin-orbit potential) are rather small. 



H. Magic Numbers for Superheavy Nuclei 



The prediction of new magic shell closures in super-heavy elements varies amongst the 



mean field models [[57 



For protons one has a competition between Z=114, 120, and 126. 
For neutrons one finds iV=172 and 184. The RMF-FR models agree in predicting a doubly 
magic 292 120 172- Precisely the same result emerges from PC-F1. This doubly magic nucleus 
is thus a common feature of relativistic models. For the density profile of 292 1 20i72, we 



observe a central depression in accordance with other mean-field approaches [ 2G , 58 , 5£ 
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we obtain, again 
deformed shell closures at Z — 104 for the 



In deformed calculations done in the way as described in Ref. 
in agreement with other relativistic models 
protons and iV = 162 for the neutrons. The nuclei in that region of the nuclear chart have 
deformations with j3 2 ~ 0.2 — 0.3. Thus also in the deformed case, these different types 
of RMF models agree well concerning their predictions of shell structure for superheavy 
elements. 
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V. QCD SCALES AND CHIRAL SYMMETRY 



QCD is widely believed to be the underlying theory of the strong interaction. However, 
a direct description of nuclear structure properties in terms of the natural degrees of freedom 
of that theory, quarks and gluons, has proven elusive. The problem is that at sufficiently 
low energy, the physical degrees of freedom of nuclei are nucleons and (intranuclear) pions. 
Nevertheless, QCD can be mapped onto the latter Hilbert space and the resulting effective 
field theory is capable in principle of providing a dynamical framework for nuclear structure 
calculations. This framework is usually called chiral perturbation theory, %PT, [O . 



Two organizing principles govern this x?T: (1) (broken) chiral symmetry (which is man- 
ifest in QCD) and (2) an expansion in powers of (Q/A), where Q is a general intranuclear 
momentum or pion mass and A is a generic QCD large-mass scale (~ 1 GeV) , which in 
a loose sense indicates the transition region between quark-gluon degrees of freedom and 
nucleon-pion degrees of freedom. Chiral symmetry is a direct consequence of the (approx- 
imate) conservation of axial vector currents. This symmetry provides a crucial constraint 
in the construction of interaction terms in the nuclear many-body Lagrangian: a general 
term has the structure ~ (Q/A) N and iV > is mandated. Higher-order constructions 
in perturbation theory (loops) will involve higher powers of {Q/A) that will, consequently, 
be smaller. This mapping from natural to effective degrees of freedom results in an infinite 
series of interaction terms whose coefficients are unknown and must be determined. 



In 1990, Weinberg |L7] introduced %PT into nuclear physics and showed that Lagrangians 
with (broken) chiral symmetry predict the suppression of N-body forces. He accomplished 
this by constructing the most general possible chiral Lagrangian involving pions and low- 
energy nucleons as an infinite series of allowed derivative and contact interaction terms and 
then using QCD energy (mass) scales and dimensional power counting to categorize the 
terms of the series according to (Q/A) N . He chose A equal to the p meson mass of 770 
MeV. This led to a systematic suppression of N-body forces, that is, two-nucleon forces are 
stronger than three-nucleon forces, which are stronger than four-nucleon forces, and so forth. 
Thus, the infinite series of interaction terms is not physically infinite. 



Following Manohar and Georgi |60| we can scale a generic Lagrangian term of the physical 
series as 
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where ip an d 7r are nucleon and pion fields, respectively, and m n are the pion decay 
constant, 92.5 MeV, and pion mass, 139.6 MeV, respectively, A = 770 MeV is the p meson 
mass as discussed above, and (<9 M , m^) signifies either a derivative or a power of the pion 
mass. Dirac matrices and isospin operators (we use t here rather than r) have been ignored. 
Chiral symmetry demands [ |6T[ 



A = / + 77-2>0 



(26) 
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such that the series contains only positive powers of (1/A). If the theory is natural |]60| , |62 
the Lagrangian should lead to dimensionless coefficients q mn of order unity. Thus, all infor- 
mation on scales ultimately resides in the Q mn . If they are natural, QCD scaling works. 

An explicit pionic degree of freedom is absent in the RMF. It has been tacitly eliminated 
in favor of an effective Hartree theory where the pion effects contribute to the various effective 
couplings. But various many-body effects are encompassed in the model parameters as well 
and may mask the underlying chiral structure. Nonetheless, it is worthwhile to classify the 
actual RMF-PC according to naturalness. Without pions, Eq. (EH) reduces to 



TfA 



f:^ 2 (27) 



and the chiral constraint Eq. (|26|) remains unchanged. 

Our test of naturalness does not care whether a particular Q mn or q„ coefficient has the 
value 0.5 or 2.0 or some other value near 1. Changing (refining) the model by adding terms 
would change all of the q mn or q„, but the same test of naturalness still applies. Adding 
new terms would simply change a specific coefficient by an amount ~ 1 (or less). That is, 
testing naturalness is largely and uniquely independent of the details, such as adding pions 
or performing more sophisticated nuclear calculations, provided the framework is given by 
Eqs. (E5f) - (|27|) while the physics is introduced via the measured observables of nuclei. 



The early RMF-PC parameterization of was tested for naturalness in [J21]J . The nine 
empirically fitted coupling constants as such span 13 orders of magnitude (ignoring dimen- 
sions). Scaling them in accordance with the QCD-based Lagrangian of |50| using Eq. (P7|), 
and taking into account the role of chiral symmetry in weakening N-body forces [|17]>H] using 



Eq. u26|), yields that six of the nine scaled coupling constants are natural. Later work f22 
refitting the model using the same Lagrangian ansatz as before resulted in two additional 
solutions where seven of the nine coupling constants are natural. These results provide 
evidence that QCD scaling and chiral symmetry apply to finite nuclei and, therefore, may 
assist in the selection of physically admissable nuclear structure interactions. However, one 



also concludes that the NHM Lagrangian |TU| may require more and/or different interaction 



terms, and this conclusion has led to our present study. It is important to note that the work 
summarized above did not test QCD, or chiral symmetry, but rather effective Lagrangians 
whose construction is constrained by QCD and chiral symmetry. 



A more extended RMF-PC adjustment was performed later [12]. This work also found 
naturalness and dimensional power counting to be extremely useful concepts in construct- 
ing realistic chiral effective Lagrangian expansions. Their expansions are based upon the 
relativistic mean-field meson models of quantum hadrodynamics (QHD) [0,0. Thus, each 
term in their Lagrangian corresponds to the leading-order expansion of one appearing in 
an appropriate |7| QHD-based meson-nucleon Lagrangian. Accordingly, their RMF-PC La- 
grangian contains nucleon densities of isoscalar-scalar, -vector, -tensor and isovector-vector, 
-tensor character, with each tensor term appearing only as a product with its corresponding 
vector term. No isovector-scalar terms appear due to their absence in the various QHD 
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approaches. In their fourth-order truncation, the best-fit set ( 16 coupling constants, uncon- 
strained search) contained fourteen natural and two unnatural coupling constants, whereas 
the worst-fit set (14 coupling constants, constrained search) is the one set containing all nat- 
ural coupling constants. Note, however, that the coupling constants of the derivative terms 
were constrained by the appropriate meson masses, as described in Sec. |11 PL Nevertheless, 



their study concludes that naturalness and dimensional power counting are compatible with 
and implied by the measured ground-state properties of finite nuclei. 

We now turn to these same considerations for the sets of coupling constants determined 
in our present study that are tabulated in Sec. III. Applying Eqs. (^6j) and (|27| ) to the 
sets of dimensioned coupling constants in Tables |TT| through |\T], and using Weinberg's 
T7| choice of the p meson mass (770 MeV) for the QCD large-mass scale A, we obtain 



the corresponding sets of QCD-scaled coupling constants listed in Table [TX|, together with 
the additional information of expansion order in A, number of coupling constants, number 
of natural coupling constants amongst them, and ratio of maximum and minimum scaled 
coupling constants in the set. The table also shows the \ 2 P er degree of freedom. The 
sets are ordered according to increasing values of this quantity. For our present work we 
require a more quantitative definition of a natural set of coupling constants than the various 
interpretations of the usual phrase "of order one" which have been applied []^,|9,21,22]: a set 



of QCD-scaled coupling constants is natural if their absolute values are distributed about 
the value 1 and the ratio of the maximum value to the minimum value is less than 10. We 



now discuss each set of QCD-scaled coupling constants appearing in Table |IX. 



A. Interaction PC Fl 

The PC-F1 interaction is the most physically realistic interaction that we have found. It 
reproduces the measured observables used to determine its coupling constants more exactly 
than any of our other interactions, as can be seen by inspection of the x\t values in Table 



IX| . Its predictive power is also better than that of the other interactions as has been 
shown in Sect. IV. The nine QCD-scaled coupling constants are all natural and the ratio 
of the maximum to the minimum is 8.92, thus satisfying our definition of a natural set of 
QCD-scaled coupling constants. So far as we are aware, this is the first complete set of 
natural QCD-scaled coupling constants, with order up to A~ 2 , that has been obtained with 
unconstrained least-squares parameter adjustment to measured ground-state observables. 



B. Interaction PC F2 

The form of the PC-F2 interaction is identical to that of PC-F1 except for the addition 
of the isovector-scalar term in Eq. (0). The most likely corresponding isovector-scalar meson 
is the S meson with a mass of 983 MeV and a relatively weak coupling constant, gg ~ 2, 
according to Machleidt [Q. Thus, its contribution is expected to be small. Nevertheless, 



the QCD-scaled coupling constant should be of order 1 if 5 meson exchange has a phys- 
ical role in the strong interaction occurring in finite nuclei in the ground state. We will 
return to this topic in our discussion of the PC-F4 interaction. Nine of the ten QCD-scaled 



24 



coupling constants of this interaction are natural whereas that of the isovector-scalar term, 
a TS = 2.34 • 10 -6 MeV -2 , is very small and unnatural, as one would expect from the above 
discussion. This small value is responsible for the relatively large ratio of 67.5 leading to 
the conclusion that this QCD-scaled set of coupling constants is not natural. This deviation 
from naturalness (here and for the following two forces) can have several reasons. There 
may be a yet undiscovered symmetry or the minimization procedure has found only a local 
minimum. 

C. Interaction PC F3 

The form of the PC-F3 interaction is also identical to that of PC-F1 except for the 
addition of the quartic isovector- vector term, Eq. (^4j) . This was done in hopes of producing 
a sign change in either of the two other isovector- vector terms, a^y or 5^y, so that their 
ratio would be positive, thus satisfying expectations based upon the 1st order expansion of 
the propagator for the p meson, as discussed in Sec. [II D| . The sign change, however, did not 
occur. Again, nine of the ten QCD-scaled coupling constants of this interaction are natural 
whereas that of the quartic isovector-vector term, 7tv — —5.4 • 10 -17 MeV -8 , is very large 
and unnatural. This large value is responsible for the very large ratio of 264.7, again leading 
to the conclusion that this QCD-scaled set of coupling constants is not natural. 

D. Interaction PC F4 

The PC-F4 interaction is built from the PC-F1 interaction by the addition of isovector- 
scalar terms that are quadratic and derivative of quadratic in the corresponding density. 
This continues the attempt with the PC-F2 interaction to address the role of the S meson 
by including both terms that are necessary to simulate the propagator. While only nine 
of the eleven QCD-scaled coupling constants are natural, and |max|/|min| is a factor ~ 10 
worse than that of the PC-F1 interaction, it is very interesting to observe that the signs of 
the two new terms are identical and thus they correctly simulate the expansion of the prop- 
agator for the S meson. Not only that, but the corresponding signs for the p meson are, for 
the first time in the present study, also identical. Thus, the expansions of the propagators 
for the two isovector mesons appearing in the PC-F4 interaction have the correct relative 
signs. Nevertheless, the maximum ratio is yet large, 103.8, leading again to the conclusion 
that this QCD-scaled set of coupling constants is not natural. We believe, however, that the 
PC-F4 interaction should be studied further. 

We conclude this section by noting that the PC-F1 interaction is one that leads to a 
physically admissable Lagrangian from the simultaneous points of view of (a) predictability 
and (b) naturalness. We have therefore demonstrated that QCD scaling and chiral symmetry 
apply to finite nuclei. 
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VI. CONCLUSIONS 



We have investigated the properties and applicability of a relativistic point-coupling 
model for nuclear structure calculations. To answer the question whether the point-coupling 
model can reach a predictive power comparable to other state-of-the-art mean field ap- 
proaches, like the RMF-FR and SHF models, we have carefully performed a \ 2 minimiza- 
tion combining two different search algorithms, and have been guided by expectations of 
naturalness in physically realistic extracted coupling constants. The resulting set of cou- 
pling constants is PC-F1 in Table |T|. It has been used to test the predictive power of the 
RMF-PC in a variety of applications ranging from saturated symmetric nuclear and neutron 
matter, binding energies in isotopic and isotonic chains to formfactor- and shell-structure 
related observables (rms charge radii, diffraction radii, surface thicknesses and spin-orbit 
splittings) and the fission barrier of 240 Pu. The net result is that the RMF-PC model with 
PC-F1 actually has reached the quality of competing approaches. In some of these compar- 
isons we discovered the influence of finite versus zero range in the models. For example, the 
density profiles of RMF-PC are generally smoother than those with RMF-FR. Like the SHF 
model, the point- coupling model naturally leads to a rather high incompressibility in nuclear 
matter, K = 260 MeV. And like the established RMF-FR forces, the point-coupling exhibits 
some unresolved isovector trends and a rather high symmetry energy in nuclear matter. The 
model performs well in deformed calculations. Also, the spin-orbit splittings are reproduced 
in a manner comparable to the finite-range models, showing that the relativistic framework 
is important here rather than the finite range. 

Attempts to extend the effective Lagrangians utilizing additional isovector terms proved 
to be elusive: the additional coupling constants can only be loosely determined with the ex- 
isting set of experimental observables. Thus the problem remains the same as in RMF-FR 
and SHF approaches, namely, that the experimental observables are very highly correlated 
with respect to the values of the coupling constants. This means that highly accurate exper- 
imental observables corresponding to large isospin are required to determine the isovector 
properties of the model more completely. 

We have been guided by naturalness in the extraction of our sets of coupling constants 
and have found that those of the set PC-F1 are all natural. In fact, so far as we are aware, 
PC-F1 is the first complete set of natural coupling constants that have been determined in 
an unconstrained search. This result, together with the predictability of PC-F1, demon- 
strates that QCD scaling and chiral symmetry apply to finite nuclei. It appears, from the 
sets PC-F2 and PC-F4, that either 5-meson exchange is not natural and is not required for 
a viable description of the strong interaction in finite nuclei, or there exists an as yet undis- 
covered symmetry. We think, however, that the PC-F4 interaction requires further study 
including possible extensions beyond eleven coupling constants (especially following new 
measurements on high-isospin nuclei) because the extracted isovector coupling constants all 
have the right signs to satisfy expectations from the expansions of their propagators. 

The point-coupling model discussed here may be viewed as a missing link between the 
established SHF and RMF-FR models. With it, one can separately investigate the in- 
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fluence of finite range versus zero range and relativistic framework versus nonrelativistic 
framework. This is important because, as we have learned, there are differences in the 
predictions from the two model classes which cannot easily be mapped onto the separate 
features of the two classes. We believe that future work should include more detailed stud- 
ies of the isovector components of the relativistic effective Lagrangians and, perhaps more 
importantly, the influence of the Fock terms via the Fierz relations. Systematic studies of 
relativistic Hartree-Fock calculations using RMF-PC will provide further linkages, on the 
one hand, with relativistic Hartree calculations using RMF-FR, and on the other hand, with 
nonrelativistic Hartree-Fock calculations using SHF. Work in these directions is in progress. 
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FIG. 1. The upper two figures show the positions of the walkers in parameter space in as 
and ay direction. For each walker, a dot is printed. The figure below on the left side shows the 
corresponding x 2 values. The other two figures indicate the size of the population and the external 
temperature. 
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FIG. 2. Errors in percent for the observables binding energy, diffraction radius, surface thick- 
ness and rms charge radius for PC-F1 (filled diamonds) and NL-Z2 (open squares) are seen on the 
left. The right panels show the absolute mean errors for the corresponding observables, where the 
dashed lines indicate the chosen relative errors AO in the fitting procedure. 



31 




I , , , , , , , 1 I , , , , , , , , , 

0.0 0.05 0.1 0.15 0.2 0.25 0.3 0.35 0.4 0.0 0.05 0.1 0.15 0.2 0.25 0.3 0.35 0.4 0.45 0.5 

Pv [fm ] p v [fm ] 

FIG. 3. The left panel shows the energy per particle versus the vector density for the forces as 
indicated. On the right side, the effective nucleon mass (upper figure) and the scalar and vector 
mean-field potentials (lower figure) are drawn as emerging from the calculations. 
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FIG. 4. Energy per particle versus neutron density for four RMF forces and the Skyrme force 
SLy6. The crosses mark data from Ref. [40| . 
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FIG. 5. Deviation in % of the calculated energies from the experimental values in spherical 
calculations of isotopic chains. Note that the scales are different for each figure. The dotted lines 
indicate the accuracy that can be demanded from the models. The experimental errors are smaller 
than the size of the symbols used in the figure. 
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FIG. 7. Deviation in % of the calculated energies from the experimental values (upper figure) 
and ground-state deformations (lower figure) in axially deformed and reflection symmetric calcu- 
lations. The errors for the binding energies are smaller than the size of the symbols used in this 
figure. The symbol with error bars indicates the measured ground-state deformation together with 



its uncertainty of i54 No [43,44|. 
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FIG. 8. Fission barrier of 240 Pu in axially- and reflection-asymmetric calculations with the 
forces as indicated. The experimental values for the ground-state deformation, the barrier height, 
and the energy of the second minimum are indicated, respectively, with an arrow, a symbol with 



error bars and three lines indicating the value and its errors. The data are taken from Ref. [ 46- 4S | . 
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FIG. 9. Surface thicknesses (lower figure), diffraction radii (middle) and rms charge radii (upper 
figure) for Sn (left) and Pb (right) isotopes emerging from spherical calculations. Note that the radii 
have been divided by A 1 ^ 3 to eliminate the liquid drop trend with mass number. The experimental 
data are from HHHI. Their 

errors are smaller than the symbols used in this figure. 
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FIG. 10. Total baryon densities of the nuclei Ca and 100 Sn emerging from spherical calcula- 
tions for the four RMF forces under consideration. 
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FIG. 11. Formfactor of the charge density of the nucleus 48 Ca for the four RMF forces under 
consideration. The experimental data, which are taken from Ref. [55], are plotted in the momen- 
tum-transfer range q = 0.35 — 3.55 [fm ], as in the original analysis. 
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FIG. 12. The percentage error in Zs-splittings for protons (left) and for neutrons (right). The 
experimental errors are smaller than the size of the symbols used in these figures. The lines serve 
to guide the eye. 
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TABLES 



observable 


error 


16Q 


88 Sr 


208p b 


Eb 


0.15 % 


+ 


+ 


+ 


rms 


0.2 % 


+ 


+ 


+ 




5.0 % 


+ 


+ 


+ 




5.0 % 


+ 


+ 


+ 



TABLE I. Observables and chosen errors AO for set 1 of nuclei used for the SA fitting 
procedure. Eb denotes the binding energy, r£ ms the rms charge radius and and ef s the proton 
and neutron spin-orbit splittings of one selected level. A + indicates an observable contributing to 
the total x 2 - For the experimental values see Ref. [ 10 1 . 



observable 


error 


03 

o o 

<© o 

1 — 1 ^ 


03 

o 

00 


% 

to 


% 

00 

in 


CO 

X 
Xj 


N 
c 

OS 


Pi 
co 
o 
o 


CO 

T 1 


co 

c 


co 


d 

CO 


Cb 
to 


3 

CO 
-< 


X 
Hh 

EN 
O 


X 

X 

c 
cs 


Ph 

CM 


Eb 


0.2 % 


+ + 


+ 


+ 


+ 


+ 


+ 


+ 


+ 


+ 


+ 


+ 


+ 


+ 




+ 


+ 




0.5 % 


+ + 


+ 




+ 


+ 


+ 




+ 


+ 


+ 










+ 




a 


1.5 % 


+ + 


+ 








+ 


















+ 




r ch 
rms 


0.5 % 


+ 


+ 


+ 


+ 


+ 


+ 




+ 




+ 








+ 


+ 


+ 


Ap 


0.05 MeV 
























+ 


+ 








A n 


0.05 MeV 
















+ 


+ 


+ 















TABLE II. Observables and chosen errors AO for set 2 of nuclei used for the Bevington fitting 
procedure. -Rdiffr. denotes the diffraction-radius, a the surface thickness, and A p and A n are the 
proton and neutron pairing gaps. A + indicates an observable contributing to the total x 2 . For 



the experimental values see Ref. [26| 
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coupling constant 


value 


dimension 


uncorr. 


error [%] 


corr. error [%] 


as 


-3.83577 • IO" 4 


MeV" 2 


2.6- 


10~ 3 


8.3 • IO" 1 


Ps 


7.68567 • lO" 11 


MeV~ 5 


2.1 • 


io- 2 


4.7 • 10° 


7s 


-2.90443 • 10~ 17 


MeV" 8 


7.2 • 


io- 2 


1.4 • 10 1 


Ss 


-4.1853 • 10~ 10 


MeV~ 4 


2.0 • 


10" 1 


2.3 • 10 1 


ay 


2.59333 • 10" 4 


MeV" 2 


3.7- 


10~ 3 


1.2 • 10° 


7V 


-3.879 • 10~ 18 


MeV~ 8 


4.4 • 


10" 1 


5.4 • 10 1 




-1.1921 • IO" 10 


MeV" 4 


6.0 • 


io- 1 


7.6 • 10 1 


axv 


3.4677 • 10~ 5 


MeV- 2 


1.2 


•10° 


1.1 • 10 1 


Sty 


-4.2 • lO" 11 


MeV- 4 


6.0 


•10 1 


1.7- 10 3 


v P 


-321 


MeV fm 3 


1.3 


•10° 


2.0 • 10° 




-308 


MeV fm 3 


1.2 


•10° 


2.3 • 10° 



TABLE III. The set PC-F1 of coupling constants resulting from the final fitting procedure. 
In columns 4 and 5 the uncorrelated and correlated errors are shown as originating from the fitting 
procedure. Note that the values for the pairing strengths have been rounded according to the error 
marigins. The total x 2 , X 2 P er point and x 2 P er degree of freedom are Xtot = 99.1, Xpt = 2.11, 
*3f = 2.75. 



coupling constant 


value 


dimension 


uncorr. 


error [%] 


corr. error [%] 


a s 


-3.835821 • 10~ 4 


MeV" 2 


1.9- 


10- 3 


8.1 • IO" 1 


/% 


7.6835 • 10~ n 


MeV- 5 


1.6- 


io- 2 


4.7 • 10° 


7s 


-2.91148 • 10~ 17 


MeV" 8 


5.2 • 


io- 2 


1.3 • 10 1 


Ss 


-4.158 • 10~ 10 


MeV" 4 


2.6 • 


io- 1 


1.8 • 10 1 


ay 


2.593511 • 10" 4 


MeV" 2 


2.9- 


10- 3 


1.3 • 10° 


7V 


-3.8234 • 10~ 18 


MeV" 8 


3.4- 


io- 1 


5.2 • 10 1 


Sy 


-1.218 • IO" 10 


MeV~ 4 


1.7 


• 10° 


6.9 • 10 1 


UTS 


2.34 • 10~ 6 


MeV" 2 


1.6 


•10 1 


2.4 • 10 3 


«TV 


3.241 • 10~ 5 


MeV" 2 


1.1 


•10° 


1.6 • 10 2 




-6.0 • 10~ n 


MeV- 4 


2.5 


•10 1 


3.7 • 10 2 


v P 


-321 


MeV fm 3 


8.7- 


IO" 1 


1.6 • 10° 




-308 


MeV fm 3 


8.1 • 


io- 1 


1.3 • 10° 



TABLE IV. The set PC-F2 of coupling constants emerging from the fitting procedure including 
the linear isovector-scalar term. In columns 4 and 5 the uncorrelated and correlated errors are 
shown as originating from the fitting procedure. The total x 2 > X 2 P er point and x 2 P er degree of 
freedom are xlt = 98-5, X 2 pt = 2-10, Xdf = 2 -80. 
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coupling constant 


value 


dimension 


uncorr. 


error [%] 


corr. error [%] 


as 


-3.835796 • 1(T 4 


MeV~ 2 


2.5- 


10~ 3 


9.9 • IO" 1 


Ps 


7.6853 • KT 11 


MeV~ 5 


2.0- 


10- 2 


5.3 • 10° 


7s 


-2.9062 • 10~ 17 


MeV~ 8 


6.9- 


10-2 


1.7- 10 1 


Ss 


-4.1797 • 10~ 10 


MeV~ 4 


2.1 • 


io- 1 


2.4 • 10 1 


ay 


2.593357 • 10~ 4 


MeV~ 2 


3.5 • 


10~ 3 


1.7-10° 


7v 


-3.8731 • 10~ 18 


MeV~ 8 


4.4 • 


10" 1 


5.9 • 10 1 


<$v 


-1.1997 • 10~ 10 


MeV~ 4 


6.9 • 


10" 1 


8.0 • 10 1 


axv 


3.549 • 10- 5 


MeV- 2 


1.2 


• 10° 


7.6 • 10° 


7TV 


-5.4 • 10~ 17 


MeV~ 8 


5.7 


•10 1 


1.8 • 10 2 




-4.0 • 10~ n 


MeV- 4 


1.1 


• 10 2 


4.0 • 10 2 


v P 


-321 


MeV fm 3 


1.4 


• 10° 


1.7 • 10° 


v N 


-308 


MeV fm 3 


1.2 


•10° 


1.3 • 10° 



TABLE V. The set PC-F3 of coupling constants emerging from the fitting procedure including 
the nonlinear term in the isovector-vector density. In columns 4 and 5 the uncorrelated and 
correlated errors are shown as originating from the fitting procedure. The total x 2 ■, X 2 P er point 



and x 2 P er degree 


of freedom are Xtot = 


98.8, x 2 P t = 2-10, xlf = 2.82. 




coupling constant 


value 


dimension 


uncorr. error [%] 


corr. error [%] 


a s 


-3.83564 • 10" 4 


MeV" 2 


2.9 • 10~ 3 


1.0 • 10° 


Ps 


7.6806 • lO" 11 


MeV" 5 


2.0 • 10~ 2 


5.7-10° 


7s 


-2.9105 • 10~ 17 


MeV" 8 


6.9 • 10~ 2 


1.9 • 10 1 


Ss 


-4.16057 • 10~ 10 


MeV" 4 


2.1 • 10" 1 


2.3 • 10 1 


ay 


2.593614 • 10~ 4 


MeV" 2 


3.5 • 10" 3 


1.5 • 10° 


7V 


-3.844 • 10~ 18 


MeV" 8 


4.4 • 10" 1 


6.8 • 10 1 


Sv 


-1.2154 • 10" 10 


MeV" 4 


6.8- 10" 1 


7.6 • 10 1 


«TS 


-5.92 • 10~ 6 


MeV" 2 


7.8 • 10° 


6.8 • 10 3 


(5ts 


-1.12 • 10~ 10 


MeV" 4 


4.1 • 10 1 


8.4 • 10 2 




3.937 • 10" 5 


MeV" 2 


1.0 • 10° 


8.9 • 10 1 




3.0 • 10~ 12 


MeV~ 4 


1.5 • 10 3 


3.1 • 10 4 


v P 


-321 


MeV fm 3 


1.4 • 10° 


1.5 • 10° 




-308 


MeV fm 3 


1.2 • 10° 


1.4 • 10° 



TABLE VI. The set PC-F4 of eleven coupling constants emerging from the fitting procedure 
including 4 isovector coupling constants. In columns 4 and 5 the uncorrelated and correlated errors 
are shown as originating from the fitting procedure. The total x' \ X 2 P er point and x 2 P er degree 
of freedom are xlt = 98.2, X 2 t = 2.09, X % = 2-89. 
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coupling constant 


value from NL-Z2 


PC-F1 


PC-F2 


PC-F3 


PC-F4 


as [Mevj 
5 S [MeV]- 4 


-1.737- 10- 9 


—o.oob • 1U 

-4.185 • 10" 10 


—o.oob • 1U 
-4.158 • 10" 10 


—o.oob • 1U 
-4.180 • 10" 10 


—o.oob • 1U 
-4.161 • 10" 10 


a v [MeV]~ 2 
5 V [MeV]- 4 


2.739 • 10~ 4 
4.502 • 1(T 10 


2.593 • 10~ 4 
-1.192 • 10" 10 


2.594 • 10~ 4 
-1.218 • 10" 10 


2.593 • 10~ 4 
-1.120 • 10" 10 


2.594 • 10~ 4 
-1.215 • 10" 10 


a T v [MeV]~ 2 
5 TY [MeV]~ 4 


3.566 • 1(T 5 
6.125 • 10" 11 


3.468 • 10~ 5 
-4.20 • 1Q- 11 


3.241 • 10~ 5 
-6.0 • lO" 11 


3.549 • 10~ 5 
-4.0 • 10~ n 


3.937 • 10~ 5 
3.0 • 10~ 12 


a T s [MeV]~ 2 
5 TS [MeV]" 4 






2.34 • 10~ 6 




-5.92 • 10" 6 
-1.12 • 10~ 10 



TABLE VII. Coupling constants from the RMF-FR interaction NL-Z2 and corresponding 



values from the RMF-PC interactions PC-F1 to PC-F4. 





PC-F1 


PC-LA 


NL-Z2 


NL3 


SLy6 


SkI3 


Po (fm 3 ) 


0.151 


0.148 


0.151 


0.148 


0.159 


0.158 


E/A (MeV) 


-16.17 


-16.126 


-16.07 


-16.24 


-15.90 


-15.96 


m* /m 


0.61 


0.575 


0.583 


0.595 


0.690 


0.577 


K (MeV) 


270 


264 


172 


272 


230 


258 


a sym (MeV) 


37.8 


37.194 


39.0 


37.4 


32.0 


34.8 



TABLE VIII. Bulk properties of nuclear matter for the forces under consideration 
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S~ \ 1 ■ /~ \ j_ j_ 

Coupling Constant 


Cl„ (PC -Fl) 


q n (PC - F2) 


ci n (PC - F3) 


c /n (PC-F4) 


Order in A 


a s 


-1.641 


-1.641 


-1.641 


-1.641 




0s 


1.443 


1.443 


1.443 


1.442 


A 1 

A 1 


75 


-2.695 


-2.701 


-2.696 


-2.700 


A —9 

A 1 


$s 


-1.061 


-1.054 


-1.060 


-1.055 


A — 2 

A z 


ay 


1.109 


1.109 


1.109 


1.109 


a n 
A u 


IV 


-0.360 


-0.355 


-0.359 


-0.357 


A -2 

A 


Sv 


-0.302 


-0.309 


-0.304 


-0.308 


A- 2 




- 


0.040 


- 


-0.101 


A 


&rs 








-1.134 


A —2 

A z 


CtTV 


0.593 


0.555 


0.607 


0.674 


A 

A u 


1TV 






-80.470 




A- 2 


5tv 


-0.422 


-0.612 


-0.404 


0.026 


A- 2 


# Ci n 


9 


10 


10 


11 




# natural 


9 


9 


9 


9 




max / min 


8.92 


67.5 


264.7 


103.8 




Xdf 


2.75 


2.80 


2.82 


2.89 





TABLE IX. QCD-Scaled Coupling Constants for Four Relativistic Point Coupling Interactions. 
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